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Abstract. The capacitance of a double- heteroj unction structure with a wide GaAs undoped 
layer embedded between two selectively doped AlGaAs barriers is calculated self-consistently 
as a function of intensity of the in-plane magnetic field. With increasing field intensity the 
capacitance initially increases and after reaching a maximum decreases toward a high field limit 
which is less than its zero field value. This behaviour is attributed to 'breathing', or charge 
redistribution, of the 2D electron gas at individual heteroj unctions due to a combination of the 
confining potential and the magnetic field. 
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1. Introduction 



Recent technological developments in the domain of semiconductors have led to tremendously 
growing interest in the (quasi) two dimensional electron gas confined to the x — y plane of 
the GaAs/AlGaAs interface by a narrow quantum well. Electrons possess only two degrees 
of freedom along the x and y direction, the motion in the z direction, perpendicularly to the 
interface, is quantized. A rough estimation of the separation between the lowest energy levels is 
about 20meV in a standard GaAs/AlGaAs heterostructure, the amplitude of the motion in the 
z direction is around lOnm. Thus the considered systems are not two-dimensional in the strict 
sense, because the wave functions always have a finite spatial extent in the z direction. 

The magnetic field B applied perpendicular to the plane of the 2D gas further quantizes the 
energy spectrum into a set of Landau levels separated by gaps with the separation Tllo determined 
by the cyclotron frequency oo = \e\B/m proportional to the magnetic field. The quantization of 
the kinetic energy of electrons by the magnetic field is responsible for characteristic magneto- 
transport phenomena, such as Shubnikov-de Haas oscillations |jl| and the quantum Hall effect 
1^]. As the effective mass m of an electron is small in GaAs, the energy quantization by experi- 
mentally available magnetic fields is of the same order of magnitude as the quantization due to 
the electrostatic potential of a quantum well. 

A magnetic field parallel to the interface does not introduce new quantum levels but strongly 
influences the electron energy structure of a quasi 2D electron system both in /c-space and in 
real space; the harmonic magnetic potential and the confining potential of a quantum well are 
combined into an effective potential for the electron motion in the z direction. Obviously, in 
such a case the system cannot be considered as a strictly two-dimensional one. 

The most pronounced effect of a parallel magnetic field is a shift of subband energy dispersion 
curves Q. The separation of levels and the density of states in branches of the energy spectrum 
increase with increasing B. As a result, it is possible to alter the individual subband population 
to the point where the highest subband is completely depopulated Q. 

It was also shown that the Fermi contours of the 2D electron gas deviate from the standard 
circular shape under the combined infiuence of the confining potential and the parallel magnetic 
field. In standard heterostructures with approximately triangular potential wells the Fermi 
contours acquire 'pear-like' shapes [^. 

The theoretical investigations of the electronic structure of a 2D system subjected to an 
in-plane magnetic field were until now restricted to simple analytically solvable models |7[ ^ , 
or to the application of perturbation theory to more realistic quantum mechanical descriptions 
of systems in a zero magnetic field . Only recently have fully self-consistent calculations of the 
energy subbands of standard GaAs/AlGaAs heterostructures in in-plane fields been published 

In this paper we focus on an interesting feature of the electronic structure of quasi 2D systems 
in parallel magnetic fields which did not attract much attention until now: the field induced 
redistribution of the electron charge density in a quantum well. This quantity is routinely 
calculated in the above mentioned self-consistent calculations but rather difficult to measure. 
New technique for making separate connections to the individual 2D electron layers of 
double-quantum wells makes it at present experimentally accessible by measuring the electric 
capacitance as a function of the in-plane magnetic field amplitude. 

In our theoretical study of the capacitance we shall consider a double-heterojunction struc- 
ture, where the wide undoped GaAs layer is embedded between two selectively doped AlGaAs 
barriers (see figure 1). With an appropriate choice of the GaAs layer width d this configura- 
tion produces two independent, non-overlapping 2D electron systems confined to GaAs/AlGaAs 
interfaces and the structure resembles a parallel plate capacitor. 

We assume that our system is symmetrical and has its centre of symmetry at z = d/2. 
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Therefore, in equilibrium, the charge distribution also is symmetrical and the electro-chemical 
potential of the left electron layer is equal to the electro-chemical potential fiji of the right 
layer. Note that in this case the electric field in the middle of the structure, at z = d/2, is equal 
to zero. 

To calculate the capacitance, we have to study the non-equilibrium state of a double-layer 
system with the left part (z < d/2) charged Ql = +Q and the right part (z > d/2) charged 
Qr = —Q- As we deal with two non-overlapping electron systems (the tunneling between 
them is completely neglected) each of them can be assumed to be in a local thermodynamical 
equilibrium, characterized by local electro-chemical potentials fj,L and hr, respectively. Then 
the capacitance C per unit area is given by 

c = M (1) 

where A^u = — 

To determine the capacitance defined by the above expression we proceed in the following 
steps. 

First, basic formulae describing the electron subbands of a quantum well subjected to the 
in-plane magnetic fields will be presented and we shall try to conclude from their shape rather 
general results independent of a detailed knowledge of the form of the confining potential. A 
simple analytically solvable model will be used to illustrate this analysis. 

Then the self-consistent calculation of the equilibrium electronic structure of the double- 
layered quantum well will be described in detail and, finally, the capacitance will be calculated, 
assuming a difference between electro-chemical potentials of individual electron layers confined 
to GaAs/AlGaAs interfaces. 



2. General formulae 

We consider a system of non-interacting electrons at each GaAs/AlGaAs interface, subjected to 
a confining electrostatic potential Vcanfi^) and to an in-plane magnetic field B = (0, By,0). The 
main difference between the electrostatic and magnetic forces is that the electric field depends 
only on the z-coordinate of an electron while the Lorentz force is also a function of the electron 
velocity components Vx and Vz- Thus, due to the presence of the in-plane magnetic field, time 
reversal symmetry is broken, and the x and z components of the electron motion couple. 

Since the translational invariance in the layer plane is preserved, the quantum mechanical 
description of the system reduces to a one-dimensional problem. In the Landau gauge, the 
momentum components Px and py are constants of motion and the wavefunction ^{r) can be 
f actor ized: 



V'(r) = ey:p{ikxX + ikyy) (fii^kA^)- 
Then the Schodinger equation for the 'out-of-plane' electron motion may be written as 



(2) 



Ei{kx) 



2m 



y^iM (^) • 



(3) 



This one-dimensional equation describes the electron motion in the effective electro-magnetic 
potential V^ff composed of the harmonic magnetic potential, corresponding to the Lorentz force, 
and the confining potential Vconf 
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Veff = {z - Zof - eVconfiz). (4) 

The centre zq of the magnetic part of the effective potential is related to the wavevector com- 
ponent kx by zq = hkx/muj. This feature of the effective potential reflects the dependence of 
the electron motion in the z direction on the x component of the in-plane electron motion. As 
a consequence, the eigenfunctions 99, ^^(5;) become /c-dependent. 

When compared with the zero field solutions, the /c-dependent eigenfunctions are modified 
by the magnetic field in two different ways. 

First, the centres of mass {z)i^k^ of ipi i.^(z) are shifted from their original positions obtained 
for B = 0. They are related to the x component of the in-plane velocity by 

{Z}i,k^ = , 5) 

mu uj 

the velocity {vx)i^kx be determined from the shape of the energy spectrum curve Ei{kx) by 

1 dEjjkx) 
^""^^^''^ = n^k^- 

Second, as the confining strengths of the magnetic potential and Vconf are added, we expect 
that the width Aj^fc^ of states (fi^k^^ {z) defined by 

Af,fc, = {z^)^^k. - {z)lk^ (7) 

will decrease with the magnitude of the magnetic field B. It can be related to the energy 
spectrum curve by equation (^) and the expression 

, 9v ^kx , , 1 dEAkx) , , 

= — ^ 8 

Both the shift of the mass centre {z)i^kx and the reduction of Aj^^^ will contribute to the 
redistribution of the electron charge density by the in-plane magnetic field. 

We illustrate the changes induced by the magnetic field in a simple analytically solvable 
model of a parabolic quantum well with the confining potential Vconf = mVL^z^ /2. In this case 
the centres of mass are shifted from the original positions {z)q = linearly with respect to kx 

I \ uj hkx 

{z)i,k^ = ^ (9) 

and the characteristic width of the wavefunction is 

\ujm J 

In the above equations uj = (u^ + Q^)^/^. Each of the equations (|9|) and ( p!o[ ) describes one of the 
above mentioned mechanisms of the charge redistribution due to the magnetic field. It follows 
from the equation (^) that the centres of mass are spread by the magnetic field until the field 
reaches a critical value 

Bc = ^- 11 

\e\ 

Then, upon futher increasing the magnetic field, the centres return to the point z = in the 
limit of infinite magnetic field. Due to the symmetry of the electrostatic potential the shifts of 
the centres corresponding to kx and —kx are symmetrical with respect to the point 2 = and 
thus the centre of mass of the whole system remains unchanged. 
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The equation (10) shows that the characteristic width of wavefunctions decreases with in- 
creasing magnetic field and in the hmit of infinite magnetic field, the wavefunctions become 
5-functions. 

One can expect this type of behaviour of eigenfunctions to be preserved for any 'reasonable' 
confining potential of a single quantum well. The centres of mass, which were localized around 
a single {z)i^kx ^ero magnetic field, occupy a finite range for B ^ 0. The occupied range first 
increases with increasing B, then decreases and finally is reduced to a point for B ^ oo. The 
positions around which {z)i^kx are localized for i? = and B ^ oo need not be the same in 
asymmetric quantum wells. 



3. Self-consistent confinement 

The standard semi-empirical model working quantitatively for the lowest conduction states of 
GaAs/AlGaAs heterostructures is used to solve the Schrodinger equation in the envelope func- 



tion approximation |12|. The envelope function is assumed to be built from host quantum states 
belonging to a single parabolic band. The effect of the effective mass mismatch is completely ne- 
glected, and the envelope functions of GaAs and AlGaAs are smoothly matched at the interface. 
Hence the Schrodinger equation has a form given by (^). 

Since we are interested in two non-overlapping double-layer electron systems confined to 
two identical GaAsAlGaAs interfaces and since the whole heterostructure is electro-neutral and 
in thermodynamical equilibrium, each of the systems can be investigated separately. We can 
restrict ourselves to a single GaAs/ AlGaAs hetero junction and demand a zero derivative of the 
confining potential Vconfiz) at the point z = d/2. 



The confining potential 



Vconf{z)=Vb{z)+Vs.c.{z) (12) 



is a sum of (i) the step function Vb{z) = Vb9{—z) corresponding to the conduction band dis- 
continuity between AlGaAs and GaAs (ii) a term describing the interaction of an electron with 
ions and (iii) the electron-electron interaction. This term should be calculated self-consistently 
and can be written as 

Vs.c.{z)=VHiz) + V^c{z). (13) 
The Hartree term Vh is determined from the Poisson equation 

f}^=\Aii^ (14) 
and for the exchange-correlation term V^c we use an expression calculated by Ruden and Dohler 



1 13 1 in a density- functional formalism 
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The conduction band offset T4 and the dielectric constant e enter our calculations as input 
parameters. (We use mks units.) 

For modulation doped GaAs/ AlGaAs heterostructures, the total charge density q{z) in equa- 
tion(p!4[) can be split into parts corresponding to concentrations of electrons, Ne{z), their parent 
donors in AlGaAs, NJ^{z), and ionized residual acceptors in GaAs, N^{z): 

e{z) = e \n,{z) - Nj{z) + N-{z)\ . (16) 



5 



The usual approximation of constant impurity concentrations is assumed. 

In our calculation we consider a GaAs/AlGaAs heterostructure with parameters N^i = 2 x 
lO^^cm"'^, Na = lO^^cm"^ and the spacer thickness w = 60nm to obtain the electron system of 
A'^e ~ 1.6 X lO^^cm"^ having one occupied subband. We took the band offset 14 = 225meV and 



the dielectric constant e = 12.9 |12|. 



As discussed in previous sections the electronic structure of a quasi 2D gas is substantially 
modified by in-plane magnetic fields both in fc-space and in real space. Figure 2 shows the 
intervals of wavevectors G (/c™*", 
fe™"^'), determined by 

Eik^n = E{krn = IX (17) 

and corresponding to occupied states (pk^ (z), together with the corresponding centres of mass. 
The {z)i^k^ curves are plotted for several magnetic fields and only for the occupied states kx S 

/t,min T^max\ 

The kx dependence of the centres is non-linear and, due to the asymmetry of the self- 
consistent confining potential, it is asymmetric with respect to the origin of coordinates. Also 
the critical values of the magnetic field and B]: for electrons belonging to the right and the 
left part of the energy spectra are different {B^ ^ 1.75T, B^ ^ 8T). 

In spite of many quantitative differences, the model with parabolic confinement gives a 
qualitative explanation of the numerical results. The electrons with kx close to A;™°^ are shifted 
into the bulk GaAs while the electrons with kx close to k!^^"' are shifted toward the interface. It 
means that these two electron systems are subjected either to a weak electric field —dVconf/dz 
in the bulk GaAs or a strong electric field —dVconf/dz at the interface. Describing the electron 
systems separately, the curve can be approximated by two straight lines (see the case of B = 4T) 
corresponding to weak parabolic confinement with $7 = and strong parabolic confinement 
with Q = Q^. Since <C the critical values of the magnetic field must fulfill B^ <^ B^ as 
follows from equation (pT|). 

With parabolic confinement the centre of mass of the whole electron system is indepen- 
dent of the magnetic field, due to the symmetry of the potential. Figure 3 shows that in the 
GaAs/AlGaAs single heterojunction the situation is different. Here the whole electron system 
is first shifted into the GaAs and then for B > B^ it is shifted back to the heterostructure 
interface. 

The distributions of electrons at i? = B^ and at B = 8T (the largest magnetic field con- 
sidered in our calculations) are compared in the figure 4(a). We can see the remarkable charge 
redistribution in bulk GaAs. At the interface, the effective electric field due to the confining 
potential Vconfiz) is so strong that the redistribution of electrons by the magnetic field is negli- 
gibly small. For this reason, the confining potential also changes remarkably only deep in GaAs 
as shown in the figure 4(b). 



4. Capacitance 

In the previous section an equilibrium state of a double- heterojunction GaAs/AlGaAs structure 
was investigated. To find its capacitance we have to study a non-equilibrium system assuming 
a difference between the local electro-chemical potentials of the left 2D electron gas, m, and 
of the right 2D electron gas, yU/?. Then the confining potential is no longer symmetrical and a 
non-zero electric field appears in the middle of the sample at z = d/2 as shown in figure 5. 

The calculation of the difference between electro-chemical potentials can be split into the 
following steps: (i) For a given value of Q the electric field in the middle of the sample, i.e. the 
derivative of the electrostatic potential at the point z = d/2, is determined by 
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d^conf ^ \e\Q .^g. 

dz y-d e 

Then the potentials V^^j{z < d/2) and V^^j{z > d/2) are calculated separately with the 
boundary condition (p!8[). To be definite, we choose the otherwise arbitrary origin of the energy 
scale in such a way that /iL is equal to zero and hr enters the self-consistent loop as a parameter 
to be determined, (ii) The electro-chemical potential fiR is obtained from the matching condition 

Kt„/W2)=Kon/W2). (19) 

In numerical practice we started with /ijj = 0, similarly as in the case of /i^, and calculate the 
difference between the resulting values V^^j and Vj^^j at z = d/2. Then obviously we have 

A/^=lKi^W2)-Kon/W2)| (20) 

and the capacitance C can be calculated from equation (|l]). We have performed the calculation 
for a number of concentrations of the non-equilibrium charge Q up to 0.5% of the equilibrium 
concentration of electrons N^, for which the value of |eA//| reaches ~ ImeV. We have found 
that the capacitance is a constant and does not vary with Q in this range of values. 

The capacitance as a function of the magnetic field is shown in figure 6. It reaches its 
maximum value at i? = and for B > B^ decreases toward a limiting value C = Coo 
for B ^ oo. The relation of the shape of this curve to the magnetic field induced charge 
redistribution is explained in figure 7. For magnetic fields B < B^ the electron systems are 
shifted closer to each other and thus the potentials V^^j and V^^j increase at the point z = d/2 
(see figure 7(a)). At the same time, since the electron concentration in the right part of the 
system is larger than in the left one, we can expect the effect of the charge redistribution to be 
stronger on V^.^j: than on Vj^^j. It means that the difference Vj^^j^ — V^^j decreases and the 
capacitance of the system increases. For B > B^ the situation is just opposite, as shown in the 
figure 7(b). 

In previous sections we mentioned that in the limit of infinite magnetic field the electron 
systems are compressed to planes close to the interfaces. It enables one to find the approximate 
limiting value of the capacitance (Coo ~ 0.57 x pF/im~^) from a model of a parallel plate capaci- 
tor with the distance between plates equal to the distance between the GaAs/AlGaAs interfaces. 



5. Summary 

The capacitance of bilayer two-dimensional electron systems as a function of in-plane mag- 
netic field has been determined by the self-consistent calculation of the electron structure of a 
GaAs/AlGaAs double-heterojunction. 

To gain better physical insight into the problem, the qualitative aspects of the behaviour 
of the system in magnetic fields have been illustrated using a simple model with parabolic 
confinement. 

The redistribution of the charge due to the combined effect of the magnetic field and of the 
shape of the confining potential was found responsible for the results obtained. 

We have also found that the self-consistent approach is unavoidable in calculation of the 
difference between the local non-equilibrium electro-chemical potentials of the two 2D systems 
which define the capacitance of a bilayer system. 

Acknowledgements- We are grateful to Prof. Allan H. MacDonald who turned our attention 
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Figure captions 



Figure 1. Schematic of the GaAs/AlGaAs double-hetero junction with non-overlaping bilayer 
2D electron systems subjected to the in-plane magnetic field B. The electron systems have one 
occupied level Eq < ji. Vq is the conduction band discontinuity, w the spacer thickness and d 
the distance between interfaces. 

Figure 2. Self-consistently calculated kx dependence of the centre of mass of the wavefunction 

in a GaAs/AlGaAs heterojunction (full lines). The discrete values of the in-plane magnetic field 
are taken from to 8T. Intervals of kx corresponds to occupied states. Dashed lines are related 
to models with parabolic confinement {Q^ <C fi^). 

Figure 3. Self-consistently calculated in-plane magnetic field dependence of the centre of mass 
of the whole system. The maximum of the curve corresponds to the critical value of the magnetic 
field ^ 1.75T. 

Figure 4. Sclf-consistcnt (a) charge distribution and (b) confining potential at B = 1.75T 
(solid line) and B = 8T (dashed line). 

Figure 5. Schematic of the GaAs/AlGaAs double-hetero junction in a non-equilibrium state 
with the left part charged +Q and the right part —Q. The local thermodynamical equilibrium 
of the bilayer electron systems is characterized by the local electro-chemical potentials /il, IJ-r- 

Figure 6. Calculated in-plane magnetic field dependence of the capacitance of the bilayer 2D 
electron systems. The maximum of the curve corresponds to the critical value of the magnetic 
field B^. 

Figure 7. Schematic of the change of the left and right part of the confining potential and 

consequently of A/i due to the magnetic fields (a) B < B^ and (b) B > B^. For B < B^, 
A// decreases with increasing the magnetic field, for B > B^ Ajj, increases with increasing the 
magnetic field. 
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